Comparing contact and dipolar interaction in a Bose-Einstein condensate 
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We have measured the relative strength Sdd of the magnetic dipole-dipole interaction compared 
to the contact interaction in a chromium Bose-Einstein condensate. We analyze the asymptotic 
velocities of expansion of a dipolar chromium BEC with different orientations of the atomic magnetic 
dipole moments. By comparing them with numerical solutions of the hydrodynamic equations for 
dipolar condensates, we are able to determine Edd = 0.159 ± 0.034 with high accuracy. Since the 
absolute strength of the dipole-dipole interaction is known exactly, the relative strength of the dipole- 
dipole interaction can be used to determine the s-wave scattering length a — 5.08 ± 1.06 ■ 10 _9 m= 
96 ±20 ao of 52 Cr. This is fully consistent with our previous measurements on the basis of Feshbach 
resonances. 
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Gaseous Bose-Einstein condensates (BEC) with 
dipole-dipole interaction (DDI) have become a fast grow- 
ing field of theoretical and experimental interest. Many 
new exciting phenomena are expected. Due to the 
anisotropic character of the DDI, most of them depend 
strongly on the symmetry of the trap. The expected 
phenomena range from modifications of the ground state 
wave function EI 01, the expansion 0,0, H^the excitation 
spectrum |a LA II, El > an d stability criteria [E 0,113 to the 
occurrence of new quantum phases in optical lattices [Tl| 
and dramatic influence on the formation of vortices and 
vortex lattices 0, 0. Dipolar BECs are now also 
discussed in the context of spinor condensates 0, 0] , 
where the combination of large spin and magnetic mo- 
ment leads to new effects like the conversion of spin into 
angular momentum |l4l Hf^ | . For all of these phenomena, 
the relative strength of the DDI compared to the contact 
interaction is a very important parameter. 
In two recent publications, we have reported on the 
generation of a BEC of chromium atoms ( 52 Cr) 
and the observation of magnetic dipole-dipole interac- 
tion (MDDI) in the BEC [h|. In the latter one, we have 
shown that depending on the orientation of the magnetic 
moments of the condensed chromium atoms with respect 
to the long axis of our optical dipole-trap, the expan- 
sion dynamics of the BEC is modified. The dynamic be- 
havior of the condensate aspect ratio after release from 
an anisotropic trap was studied experimentally and com- 
pared to numerical calculations based on the description 
of the chromium BEC by superfluid hydrodynamic the- 
ory including dipole-dipole interaction [5J, |8|, Il9j . The 
observed behavior showed an excellent qualitative agree- 
ment with the theoretical prediction. In this paper, we 
discuss a method that allows to determine the strength 
of the MDDI compared to the contact interaction with 
an accuracy on the percent level. 



Since the absolute strength of the MDDI is known, 
one can use a measurement of the relative strength of 
the MDDI to determine the s-wave interaction strength, 
which is proportional to the s-wave scattering length 
a. Many different techniques have been used to deter- 
mine the s-wave scattering lengths of ultra-cold atoms 
but most of them come along with large error bars 
- often because the number of atoms enters the mea- 
surement. Examples are the 23 Na scattering length in 
| F = 1, m F = -1 ) of a Na = 92 ± 25 (i.e. 27% error) 
determined from thermalization measurements |2(il |. and 
aNa = 65 ± 30 (i.e. 46% error) from the measurement 
of the mean field energy of a BEC |2l|. The scattering 
length of metastable He* anc* = 16 ± 8nm was deter- 
mined from the mean field energy by analyzing the size 
of the BEC where the error of 50% stems from an un- 
certainty of the number of atoms in the condensate. Our 
first determination of the chromium scattering length j2^| 
was based on cross-dimensional thermalization measure- 
ments [24j and resulted in act =170±39ao. The error was 
mostly due to an uncertainty in the density and atom- 
number determination. We will show in this paper, that 
a measurement of the relative strength of the magnetic 
dipole-dipole interaction in a BEC can be used to obtain 
precise values for a without such a strong dependence on 
the determination of the number of atoms. 
The interaction energy of two magnetic dipoles separated 
by the distance f is given by 
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where the strength of the dipole-dipole interaction is 
measured by the pre-factor of Udd and the orientation 
of the dipoles is parallel to an external magnetic field 
B. This strength can be compared to the coupling con- 
stant g of the s-wave interaction 
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and is measured by the dimensionless dipole-dipole 
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strength parameter 



gives rise to the mean-field potential Q, Q 
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It is chosen such that a homogeneous condensate is un- 
stable if Sdd > 1 in a static magnetic field |25j . 
In contrast to the s-wave interaction which can be under- 
stood as a local, contact-like interaction the dipole- 
dipole interaction is long-range and anisotropic. In a 
condensate with density distribution n{f) = \4>{f )| 2 , it 



*dd(f)= I U dd (r-r'Mr')\ 2 d 3 r' 



(4) 



The integral in (0} reflects the non-local character of the 
interaction. If this interaction in addition to the contact 
interaction is taken into account, the well known Gross- 
Pitaevskii equation gets the form 



' — V 2 + U ext (r)+g\<j>(r,t)\ 2 + / U dd (f - f ')\cb(r ', t)\ 2 d\ ' 0(r, t) 
2m ' 



(5) 



O'Dell et al. have shown in [8j, that even under the influ- 
ence of the dipole-dipole mean field potential 3>dd(f), the 
density distribution has the shape of an inverted parabola 
in the Thomas- Fermi limit. Like in the case of pure con- 
tact interaction, a wave function of the form 
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is a self consistent solution of the superfluid hydrody- 
namic equations j2(J derived from the Gross-Pitaevskii 
equation ©, even in presence of dipole-dipole interac- 
tion |S, 19]. R x , Ry, and R z are the Thomas- Fermi radii 
of the condensate. The anisotropy of the dipole-dipole in- 
teraction manifests itself in a modification of the aspect 
ratio of the trapped condensate |l9l |25| . This anisotropy 
also reveals during the expansion of a dipolar conden- 
sate nm. 

In the following we will determine the dipole-dipole 
strength parameter e dd by analyzing the dynamic behav- 
ior of the Thomas-Fermi radii R(t) of expanding dipo- 
lar condensates. The experimental apparatus and tech- 
niques that are used are described in detail in H, [27]] • By 
applying a small homogeneous external field (~11.5G), 
oriented either along the y— or z— axis, shortly (~7ms) 
before releasing the condensate from the trap, we ob- 
tained two sets of measured radii of a ballistically expand- 
ing condensate with different alignment of the atomic 
magnetic moments. The trap from which the conden- 
sate was released was a crossed optical dipolc trap that 
was elongated in z-direction with trap parameters of 
w x = 2tt 942 Hz, uj v = 2tt 712 Hz, and u z = 2n 128 Hz. 
In the asymptotic limit of long times of flight, which is 
governed by a collisionless and potential free (except for 
gravity) ballistic flight, the radii of the cloud can be pa- 
rameterized as 



Ri(t)=R*+v*t, 



(7) 



where the index i — [x,y,z] indicates the direction of 
expansion that is considered. It is worthwhile to mention 



polarization 


v*[10- 3 m/s] 


C [m 4 / 5 ] 


no dipoles 


8.528 


0.0488 


y— polarization 


9.085 


0.0519 


2— polarization 


8.283 


0.0474 



TABLE I: Asymptotic velocity in y— direction and corre- 
sponding proportionality constant C calculated numerically 
for the case of vanishing dipole-dipole interaction Edd = 0, 
and for Edd = 0.148 and polarization along y and z. Veloci- 
ties calculated for 30000 atoms, and a = 103ao. 



that the initial values Ri(Q) — R* are not the Thomas- 
Fermi radii Ri. Note that the radii Ri(t) as well as the 
asymptotic velocities v* of the expansion for long times 
(t 3> 1/lu) depend on the direction in which the atoms 
are polarized. As shown in , they are thus proportional 
to (iVa) 1 / 5 . In particular, 
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with a constant of proportionality C that only depends 
on the known or measured quantities that determine the 
chemical potential, i.e. the trap parameters uj x ,LUy,uj z 
the atomic mass m and a small contribution of e dd . Us- 
ing the hydrodynamic theory of an expanding dipolar 
condensate , the asymptotic velocity for a certain num- 
ber of atoms and scattering length can be easily calcu- 
lated numerically. Table [I] shows the expected asymp- 
totic velocities v* = C ■ (30000 • 103a Q y/ 5 in y-direction 
and the corresponding values for C, calculated for 
pure contact interaction (e dd = 0), y— polarization and 
z— polarization. The numbers are calculated for the mea- 
sured trap parameters, 30000 atoms, and a scattering 
length of a = 103ao 28]. The scattering length of 
103ao corresponds to a dipole-dipole strength parame- 
ter of e dd = 0.148. 

To determine the asymptotic velocity, we use the con- 
densate radii R(t) measured in a time-of-flight series. 
Since the number of camera pixels that are covered by the 
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FIG. 1: Measured dependence of the condensate radius R y (t) 
on the time of flight, (a) dipoles aligned along the y— axis 
(B\\y), average number of atoms was 29000 (b) dipoles aligned 
along the z— axis (-Bp), average number of atoms was 31000. 
Open circles: measured radii, crosses: measured radii re- 
scaled using (0; solid black line: linear fit to the re-scaled 
radii. 



condensate is much larger in the direction of fast expan- 
sion (y— direction in our setup), one can expect the most 
accurate results when considering R y (t). Figure ^ shows 
the dependence of the condensate radius R y (t) with po- 
larization along y (left figure) and z (right figure) on the 
time of flight for 67 different expansion times. Because 
the number of atoms is fluctuating during such a series 
of experiments, the radii fluctuate due to their depen- 
dence on To get rid of these fluctuations, we divide 
each measured radius R y by the fifth root of the number 
of condensed atoms N in the corresponding experiment 
and multiply them with the mean value of the fifth root 
of the atom-numbers < TV 1 / 5 > in all experiments. 
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fn this way we get a series of time dependent radii which 
are now independent of the atom-number. Open cir- 
cles in Fig. n represent the measured R y {t), crosses with 
error-bars mark the re-scaled R(t) which show much less 
fluctuations. A linear fit to the re-scaled data for times 
larger than 3 ms to focus only on the asymptotic behav- 
ior, yields v* — 9.56 ± 0.24m/s for y— polarization. For 
z— polarization, we get v* — 8.78 ± 0.12m/s. 
By using the above re-scaling, the errors Av* — 
±0.24m/s and Av* = ±0.12m/s in the fitted slope v* 
for y— and z— polarization, respectively, do not contain 
fluctuations of the atom-number anymore. 
If we consider only the case of y— polarization, invert 
equation (JHJ) and insert the fitted velocity v* and the 
constant C y from table |U we get 



(10) 



for the scattering length. The error on this measurement 
consists of two contributions: first the fitted asymptotic 
velocity comes with an error Av* due to such kind of 
noise on the data that is not correlated with the number 
of atoms. Since v* appears in the fifth power in a, Av* 
appears with a factor of 5 in Aa. Second, the mean value 



of the number of atoms has an uncertainty, mainly due to 
an uncertainty in the detuning of the probe beam. Since 
a detuning from resonance can only lead to an underesti- 
mation of the number of atoms, the error in the scatter- 
ing length a caused by this uncertainty is only towards 
smaller values of a. We estimate a maximum detuning of 
A5 pro b e = ±0.25r which leads to an estimated error in 



the number of atoms of AN/N 
error in a is then 



-0.25%. The relative 



Aa A < N > Av* , s 
= + 5 = -0.25 ±0.075. 11 

a < N > v* 

Hence the scattering length of 52 Cr determined with this 
method is 



a Cr = (l38^jj)a . 



(12) 



For z-polarization, we get a consistent value of eto = 



(13312> 



-4sj a o- Due to the relatively large systematic error 

in the number of atoms, this way to determine the scat- 
tering length yields only quite inaccurate values, typical 
for condensate expansion experiments. In the following, 
we will use the full set of data from both polarisations to 
determine the scattering length with much higher accu- 
racy and independent of the number of condensed atoms. 
We use the two rescaled asymptotic velocities 
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< N 1 / 5 > 



v*(B\\y) and v*(B\\z) for polarization along y and z, re- 
spectively to determine Edd by analyzing their ratio. To 
first order in Sdd (in the expected range of Edd, higher 
orders are negligible), the ratio has the form 



V*y (B\\Z) 



= 1 + De dd - 



(14) 



It depends only on the asymmetry introduced by the 
dipole-dipole interaction because the contribution of the 
s-wave scattering to the total energy is independent of 
the polarization. D is again a numerical constant. If we 
use the measured asymptotic velocities, we obtain 



e dd = 0.159 ±0.034 



(15) 
0.148 



in very good agreement with the value of Edd 
that one would expect for acv = 103ao- 
In turn, since the dipole-dipole interaction strength can 
be exactly calculated from J3J , this result can be used to 
determine the scattering length: 



12irh 2 t 



= 96 ±20 o . 



(Id 



This result is in excellent agreement with the value of 
103±13a that has been obtained by comparing the mea- 
sured positions of Feshbach resonances in chromium col- 
lisions with multichannel calculations 28]. Furthermore, 
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method 


year 


scattering length acr [ao] 


relative error 


cross-dimensional thermalization [23 


2003 


170 ± 39 


23% 


Feshbach resonances \28j 


2005 


103 ± 13 


13% 


condensate expansion in one direction (this 


paper) 2006 


138 and 133 +4° a 


+8 Of 
-33 /0 


dipolar expansion (this paper) 


2006 


96 ±20 


13% 



TABLE II: Comparison of experimental values of the 52 CV s-wave scattering length. 



the relative error of ±20% makes this way to determine 
the scattering length a comparably accurate method. 
Table [H] shows a comparison of the values of the s-wave 
scattering length of chromium measured with different 
methods. 

In conclusion, we have measured the relative strength 
of the magnetic dipole-dipole interaction in a chromium 
BEC. The relative strength parameter was measured 
£dd = 0.159 ± 0.034 by analyzing the expansion of a 
chromium condensate with different polarization after re- 
lease from an anisotropic trapping potential. This re- 
sult was used to determine the s-wave scattering length 
of 52 Cr a = 96 ± 20 ao in excellent agreement with 
the results of theoretical analysis of measured Feshbach 
resonances (a52 Cr = 103 ± 13ao) [2^. In contrast to 
many other methods that are commonly used to deter- 
mine the s-wave scattering length, this method does not 
depend on the accuracy of the atom-number determi- 
nation. Furthermore, unlike other methods that deliver 
results with the same or even better accuracy like Fesh- 



bach resonance measurements or photoassociation spec- 
troscopy [2^, H3, 0, it does not require knowledge of 
any details of the molecular potentials. We expect it to 
be well suited to determine the scattering length close to 
a Feshbach resonance with high accuracy, especially for 
small scattering lengths where the dipole-dipole interac- 
tion becomes as important as the contact interaction. 
The excellent agreement between experimental results 
and theory constitutes a confirmation of the theoretical 
approach that is used to describe the dipolar BEC. 
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